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ABSTRACT 

Well over 700 exoplanets have been detected to date. Only a handful of these have been observed directly. 
Direct observation is extremely challenging due to the small separation and very large contrast involved. Imaging 
polarimetry offers a way to decrease the contrast between the unpolarized starlight and the light that has become 
linearly polarized after scattering by circumstellar material. This material can be the dust and debris found in 
circumstellar disks, but also the atmosphere or surface of an exoplanet. 

We present the design, calibration approach, polarimetric performance and sample observation results of the 
Extreme Polarimeter, an imaging polarimeter for the study of circumstellar environments in scattered light at 
visible wavelengths. 

The polarimeter uses the beam-exchange technique, in which the two orthogonal polarization states are 
imaged simultaneously and a polarization modulator swaps the polarization states of the two beams before the 
next image is taken. The instrument currently operates without the aid of Adaptive Optics. To reduce the effects 
of atmospheric seeing on the polarimetry, the images are taken at a frame rate of 35 fps, and large numbers of 
frames are combined to obtain the polarization images. 

Four successful observing runs have been performed using this instrument at the 4.2 m William Herschel 
Telescope on La Palma, targeting young stars with protoplanetary disks as well as evolved stars surrounded by 
dusty envelopes. In terms of fractional polarization, the instrument sensitivity is better than 10~^. The contrast 
achieved between the central star and the circumstellar source is of the order 10~^. We show that our calibration 
approach yields absolute polarization errors below 1%. 

Keywords: instrumcntation-polarimetry - circumstellar disks - exoplanets 

1. INTRODUCTION 

More than 700 exoplanets have been detected to date [l] . Almost all of these have been detected using indirect 
techniques such as radial velocity surveys and transit detections. The information that can be gathered about the 
planets from these detections is however limited. For the majority of the planets, detected through radial velocity 
measurements, only the orbital eccentricity and distance and a lower limit of the planetary mass can be obtained. 
Transits yield the planet radius, mass, orbital inclination and in a few cases the upper atmospheric composition 
[2j|6|. Moreover, these two detection methods are subject to selection effects: radial velocity detections favor 
high-mass planets in close orbits. Transit detections are also biased towards short-period orbits and can only 
detect planets in orbital planes intersecting our line of sight to the star. 

To further characterize exoplanets we must be able to directly detect light from them and separate it from 
the light of the central star. Direct imaging of exoplanets is however a very challenging prospect. Not only is the 
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Table 1. ExPo key parameters. 



Wavelenght range (nm): 500 - 900 

Field of view (" per beam): 20x30 

Image scale {" /pix): 0.078 

Detector type: Andor iXon EMCCD 

Size (pixels): 512x512 

Pixel pitch (/im): 16 



separation between a star and its companion very small; the ratio of their intensities, the contrast, is extremely 
high. As a consequence, only a handful of exoplanets have so far been imaged directly (e.g. [Tjjs]). These 
direct detections have all been performed in the (near) infrared, imaging warm, self-luminous planets around 
young stars. To observe and characterize cooler planets or even contemplate imaging of rocky planets, we require 
observations of starlight at visible wavelengths reflected by the planetary atmosphere or surface. Unpolarized 
light from the central star that is reflected by exoplanets becomes linearly polarized by the reflection process (e.g. 
[9]). Imaging polarimetry is therefore an ideal way to discriminate between the polarized light from circumstellar 
environments and the unpolarized light from the nearby central star. The use of polarimetry to detect and study 



extra-solar planets has already been investigated by [10 ■ 12 



We have built the Extreme Polarimeter (ExPo) with the aim of advancing this polarimetric technique. This 
imaging polarimeter has been developed with the goal of achieving an imaging polarimetry contrast better than 
10~^ unaided by Adaptive Optics. The first concepts for this instrument have been presented by jlSi, the initial 



design by 14 . Three successful observing runs have been performed using this instrument so far, all at the 
Nasmyth visiting focus at the 4.2 m William Herschel Telescope (WHT) on La Palma. The targets observed 
have included stars with circumstellar disks and protoplanetary disks or debris disks. In this paper we present 
the instrument design, polarimetric performance and sample results. 

We present the measurement principle, beam-exchange polarimetry, in section 2. Section 3 describes the 
instrument design and section 4 the polarimetric calibration approach. Results from the first observing runs 
illustrating the instrument performance and scientific potential are presented in section 5. Our conclusions are 
presented in section 6. 

2. BEAM-EXCHANGE POLARIMETRY 

Polarization measurements are inherently differential: In the case of linear polarization, two measurements taken 
through orthogonal polarization filters are subtracted. An unpolarized source will be equally present in both 
measurements and vanish in the difference, while a polarized source will remain. It is necessary to make the 
two measurements simultaneously, as temporal changes caused by atmospheric seeing will otherwise appear as 
false polarization signals. A polarizing beamsplitter can be used to do this. However, differential aberrations 
between the two beams will still appear as false signals. To solve this, a polarization modulator is used to 
alternately rotate the polarization by 90°, effectively swapping the polarization information that is passed by 
the two beams. This so-called beam-exchange technique has been developed for sensitive solar polarimetry in 



the presence of atmospheric seeing 15 . 16 has shown this technique to limit errors to 2nd-order effects and to 
be able to achieve polarization sensitivities better than lxl0~^ of the total intensity. The technique is depicted 
schematically in Figure 1. 

To reduce the effects of atmospheric seeing, we take long sequences of many short-exposure frames, combining 
the frames in sets of 4 and averaging over the sets. The short exposure times (typically 28 msec) ensures that 
the instantaneous PSF is relatively unchanged between the frames that are subtracted which would otherwise 
lead to false polarization signals. The alignment of frames within a set is performed based on the PSF, and may 
change slightly between one set and the next. Using a polarization modulator has a further advantage: The 
filtering of the polarization state to which the instrument is sensitive takes place at the modulation element, 
making the instrument insensitive to instrument-induced polarization downstream of the modulator. It is thus 
desirable to place the modulator as far upstream in the instrument as possible. 
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Figure 1. Diagram demonstrating the dual beam-exchange technique being used to detect an exoplanet. The modulator 
alternately aligns the polarization orientation of interest with two orthogonal directions that are then separated by the 
polarizing beamsplitter. Both beams are imaged simultaneously by a camera. Subtracting the resulting images cancels 
out the large unpolarized signal while the (partially) polarized planet signal remains. 

3. INSTRUMENT DESIGN 

The ExPo design was created to perform five distinct but linked functions, each of which will be discussed in 
detail below: 

Polarimetry: The core of the polarimeter consists of the polarization modulator and the analyzer. For dual- 
beam polarimetry, the analyzer needs to be a polarization beamsplitter. We have chosen a cube beamsplitter 
over crystal-based alternatives such as a WoUaston or (modified) Savart plate. The advantage of the cube is that 
the polarization beamsplitting is achromatic and does not suffer from crystal-induced astigmatism. The price we 
pay for this is the lower extinction ratio. Moreover, the extinction ratios is not the same for the reflected and the 
transmitted beams. ExPo currently uses a standard polarizing beamsplitter cube offered by Halle that provides 
an extinction ratios of 1:2000 and 1:40 for the transmitted/reflected beam over the 600-900 nm wavelength range. 

The choice of polarization modulator has been driven by the fact that we wish to observe using short exposures 
in order to minimize the seeing variation during individual exposures and thus need to switch polarization states 
rapidly between exposures to achieve a high duty cycle. This means a mechanically rotating half-wave plate is 
not an option. We have chosen to use a Ferro-electric Liquid Crystal (FLC) device, which offers state transition 
times below 100 /isec. This is significantly faster than twisted nematic cells, for which the transition time is 
several milliseconds. Another difference is that with FLCs, the orientation of the fast axis is switched, while 
in nematics, it is the retardance that is modulated. A problem with FLCs (as with the nematics) is that by 
definition, the device is a half-wave retarder at a specific wavelength, and so is chromatic. It is possible to combine 
several waveplates in an achromatic combination [17' and we intend to do this using three FLCs. However, other 
attempts to do this have not yielded encouraging results and so far we have used a single FLC centered on 
750nm. 

The orientation of the FLC with respect to the polarizing beamsplitter determines the orthogonal polarization 
components (e.g. Stokes ±Q) the instrument is most sensitive to, namely the the polarization angles that receive 
the maximum modulation. Polarization components at 45° to these angles (e.g. Stokes ±U) will have minimal 
modulation, and thus the instrument will be insensitive to them. The FLC is therefor mounted on a rotation 
stage so the full linear polarization field can be measured. 

ExPo has been designed to operate from a gravity-invariant position in the Nasmyth focus of the William 
Herschel Telescope. At this position, a significant amount of linear polarization is introduced by the 45° M3 
reflection. The ExPo design includes a tilted glass plate to compensate this unwanted polarization. A control 
loop adjusts the tilt continuously after measuring the introduced polarization by comparing subsequent images 
taken at opposite modulator. The tilt assembly is mounted on a rotation stage that keeps the tilt axis aligned 
with the M3 orientation as the telescope elevation changes. In order to perform polarimetric calibration, a linear 
polarizer can be inserted into the beam and rotated to any desired angle. 
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Coronagraphy: To reach high polarimetric sensitivity with a given dynamic range of the detector, it is desirable 
to suppress the peak of the steUar PSF. An added benefit is that the stray hght in the instrument, and the 
possibility of it becoming polarized, is reduced. To this end, the ExPo design incorporates a coronagraph. We 
have opted for a classical Lyot-coronagraph consisting of a mask in the focal plane and a Lyot stop placed in 
the pupil. A large number of other coronagraph concepts exist but most tend to either introduce chromatic 
effects or be incompatible with polarimetry. Using a focal mask etched on a glass substrate allows us to use 
semi-transparent masks. The transmitted stellar peak can be used for precise guiding. The ExPo focal mask can 
be translated to place coronagraph dots of different sizes in front of the star. The Lyot stop includes masking of 
diffraction by the telescope spiders, and is mounted on a rotator stage that keeps the mask aligned with changing 
telescope elevation. 

Imaging: In order to record both beams simultaneously and avoid synchronization issues associated with 
multiple detectors, ExPo relies on a single detector onto which both beams are projected. Two right-angle 
prisms mounted onto the polarizing beamsplitter are used to achieve this, as shown in Figure 2. 




Figure 2. After the beamsplitter, the beams are projected onto the detector by slightly rotating the two right-angle 
prisms. 

The detector itself is rotated by 45° to align with the symmetry axis created by the beamsplitting cube. 
A disadvantage of this design is the slight projection angle introduced in the images. This effect is however 
symmetric; The images can still be combined. As the instrument operates without the use of an Adaptive Optics 
(AO) system to correct for the atmospheric seeing, we wish to take large numbers of short exposures, with 
exposure times of the order of the speckle lifetime. Yet we also wish not to be limited to observing extremely 
bright targets by the detector readout noise. This has led to the selection of a detector using Electron-Multiplying 
CCD technology, where a signal can be amplified on the detector before being read out. The detector is an Andor 
iXon with a resolution of 512 x 512 pixels and a maximum framerate of 35 Hz. The design of the imaging optics 
also includes two fairly slow (//# 10) lenses, with most optical elements placed in the collimated beam between 
them. The focal ratio demagnifies slightly to create a field of view of 20 x 30 arcseconds on each half of the 
detector. A smaller field of view can be obtained by increasing the focal length of the second lens. This lens is 
placed directly before the polarizing beamsplitter. 

Atmospheric dispersion compensation: The dispersion introduced by the atmosphere needs to be com- 
pensated before the coronagraph mask placed in the prime focus. A longitudinal ADC design has been chosen, 
primarily because any polarization introduced will be stationary in the telescope frame of reference. As with the 
polarization introduced by the Nasmyth flat, this can be compensated by the polarization compensator. The 
design consists of two wedge prisms, one fixed and one that can translate along the optical axis to increase the 
distance between the prisms, increasing the dispersion angle compensated. The complete ADC is mounted to 
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the telescope derotator flange that can rotate to keep the prisms ahgned with the plane of dispersion as the 
telescope elevation changes. For this design it was possible to use two identical off-the-shelf 3-inch prisms with 
a wedge angle of 3°. 

Wavelenght selection: A filter wheel with 12 50 mm circular filter positions is used to perform ExPo obser- 
vations with an array of broad- and narrowband filters. 

Figure 3 presents the instrument block diagram below showing the integration of these five design elements. 



A" 



J" 




Detector 



Figure 3. Schematic diagram showing the different ExPo components. The occulting mask in the prime focus and the 
aperture mask in the pupil image together form the coronagraph that suppresses the light from the star. The core of 
the polarimeter consists of the Ferro-electric Liquid Crystal (FLC) acting as the modulating element and the polarizing 
beamsplitter that splits the beam into its orthogonal polarization components. 

In principle it is desirable to place the polarization modulator as far upstream in the instrument as possible. 
Any instrumental polarization introduced downstream will not be modulated and so cancel out in the differential 
imaging. The modulator however needs to be placed after the polarization compensator. As the compensators 
tilted glass plate introduces a slight and variable beam shift, it must be placed after the pupil mask, which 
requires precise alignment. This is the reason the modulator ends up relatively far downstream in the design. 

Mechanically, the design is intended to be flexible, with components mounted and aligned individually on 
the optical table at the Nasmyth visiting instrument focus of the William Herschel Telescope. Much of the 
mounting hardware is off-the-shelf, with custom mounting hardware produced by our in-house workshop. After 
each observing run, the instrument is shipped back to our institute for testing and improvements in our optical 
lab. Figure 4 shows the instrument as mounted at the telescope. 

4. DATA REDUCTION AND POLARIMETRIC CALIBRATION 

The raw product of ExPo measurements are image series in which the even and uneven frames correspond to the 
two different states of the FLC polarization modulator. Each raw image has a left and right sub-image formed 
by the two beams of the polarizing beamsplitter. 
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Figure 4. Image showing the instrument as installed on the WHT Nasmyth port optical table during one of the observing 
runs. 




Figure 5. Raw ExPo frame showing the double image created by the two beams. 

The FLC itself can be rotated to change the polarization direction to which the instrument is sensitive. The 
angle at which the instrument is sensitive to ±Q (in Stokes notation) is offset 22.5° from the angle for zLU. 
For redundancy, ExPo observations are performed at four different FLC angles: 0°, 22.5°, 45°, 67.5°. With the 
combination of two FLC states (A, B) and the left/right sub- images, the beam-exchange technique thus yields 
four different types of images at every FLC angle. We denote these intensity images by: 
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Table 2. ExPo Photon budget (650 nm) 



Element 


Transmission 


Atmosphere (V-band) 


0.90 


Central obscuration 


0.94 


Telescope primary 


0.90 


Telescope secondary 


0.90 


Nasmyth flat 


0.90 


ADC 


0.98 


Focal mask (no dot) 


0.95 


CoUimating lens 


0.98 


Filter (Schott OG590 glass) 


0.90 


Lyot stop 


0.75 


Polarization compensator 


0.93 


Polarization modulator (PLC) 


0.85 




98 


Beamsplitter (reflected) 


0.48 


Beamsplitter (transmitted) 


0.42 


Right-angle prisms 


0.98 


CCD Quantum Eflaciency 


0.93 


Total (reflected beam) 


0.13 


Total (transmitted beam) 


0.11 


Photons, TOy = star: 




Entering 4.2 m telescope: 


1.12x10" ph/s 


Detected (reflected beam): 


1.4x1010 ph/s 


Detected (transmitted beam): 


l.SxlQio ph/s 



For a general polarimeter, a minimum of four independent signals must be measured in order to determine 
the four elements of the Stokes vector. The signal matrix X relates the measured signals S to the input stokes 
vector I: 

S = X-I (1) 

If X can be determined from a calibration source, its inverse can be used to obtain the Stokes elements for 
the observed targets: 

I = Y-S with Y = X-i (2) 

ExPo does not measure circular polarization, so we are only interested in the first three elements of the Stokes 
vector: /, Q and U. For this, four measurements are more than sufiicient. The two beams, two FLC states 
and four FLC angles combined provide 16 difi^erent measurements from which four must be chosen. However, 
choosing the A/B and left/right combination of images for a single FLC angle does not work. As the FLC 
exchanges the beams between the two modulation states, the A-left images will to a large degree contain the 
same polarization information as the B-right images. The same goes for the B-left and A-right pair, resulting 
in an ill-conditioned signal matrix X. Another option is to combine all four FLC angles for a single choice of 
FLC state and beam, e.g. A-lcft. These four measurements are sufficiently independent and matrix X will be 
well-conditioned and invertable. But this approach results in another problem: As the observation at a single 
FLC angle may take up to an hour, observational conditions will have changed too much between the first and 
the last FLC angle for them to be combined in a meaningful way. Furthermore, as the telescope tracks the target 
across the sky, unpolarized artifacts such as ghosts may change position in the image and will then no longer 
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Table 3. Dome flat polarization after each calibration iteration and the effect of the Polarization Compensator (PoCo). 



Calib. 


PoCo Off 


PoCo On 


iteration 


Q (%) 


U(%) 


Q(%) 


U(%) 


1 


0.0 


0.0 


0.0 


0.0 


2 


-0.02765 


3.22021 


0.00437 


0.00243 


3 


-0.02740 


3.21953 


0.00437 


0.00243 


4 


-0.02740 


3.21953 


0.00437 


0.00243 



vanish from the result. To overcome these issues, we start by combining the images for a single FLC angle in 
the so-called double-difference: 

As the left and right sub-images are located on the same raw image, it is essential that they are carefully co- 



aligned before the above subtraction takes place. 18 discusses in detail the processing of raw images performed 



to generate the double-difference images. Unpolarized signals drop out in this measurement, for which 16 
has shown that all first-order effects drop out. A reduced signal matrix now relates the two double-difference 
measurements for FLC angles 22.5° degrees apart to the linear polarization elements of the Stokes vector: 

^ ^Stl ] =X. ( ) (4) 



\ Jno j \ 

The parameters of X are determined by taking dome flats with the calibration polarizer at different angles. 
A fitting routine is used to fit the response parameters a, b for each double-difference signal: 

= ( « O • ( ^ ) (5) 

The initial fit assumes the light entering the calibration polarizer is unpolarized. In reality this is not the case: 
The dome lights may be polarized and the Nasmyth flat will certainly introduce polarization. We iterate the 
fitting by using the first fit to calibrate dome flats taken without the calibration polarizer and estimate the real 
input Stokes vector. A new fit is then performed using this estimate. Errors are estimated by back-calibrating 
the flats taken with the calibration polarizer. Two iterations are sufficient to reduce the errors to a few percent. 
For real observations, the instrumental polarization introduced upstream of the calibration polarizer is removed 
by the polarization compensator. Table 3 shows the increase in accuracy of the consecutive calibration iterations 
and the effect of the polarization compensator. 

The current ExPo setup is quite chromatic, using only a single FLC and not an achromatic, Pancharatnam 
configuration. As we use dome-flats with a specific spectrum for our broad-band calibration, the calibration 
may not be correct for targets that have a substantially different spectrum. This is certainly the case for young 
stars peaking towards the blue, such as AB Aurigae, as the dome flat lamps are fairly red iTefi = 3000K). 
To overcome this problem, we have created a model of ExPo including the spectral response of all wavelenght- 
sensitive components, in particular the FLC, the polarization beamsplitter and the QE-curve of the detector. 
Using this model, we are able to generate synthetic calibration flats for any desired spectrum. The model has 
been calibrated to so that integrated dome flats generated with the dome flat spectrum are within 1% of the real 
flats. 

As unpolarized signals drop out in the double difference. Stokes / cannot be determined from the double 
difference, but is instead found from: 

^ — 2^^A,Uft + ^A.right + ^B,l^ft + '^'fj^^ighJ (6) 
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Note that this yields the cahbrated intensity in polarimetry terms but does not take attenuation by the 
atmosphere and optical transmission losses into account. A photometric intensity calibration must be performed 
using standard stars. The final step in the data reduction is the correction for sky rotation during the long 
observation sequences. This is achieved by first obtaining the double difference measurements from subsets 
corresponding to 30 seconds of observations. Each block of frames is physically rotated and calibrated while 
correcting for the rotation of polarization frame of reference. These blocks are then added together to obtain 
the final images in Stokes /, Q and U. 



5. INSTRUMENT PERFORMANCE 

Three successful observation campaigns have been conducted using the ExPo instrument at the William Herschel 
Telescope on La Palma. A short fourth run was conducted with a smaller field of view of 10" x 15" , but this 
did not produce good results, partly due to adverse meteorological conditions. The observations have initially 
focused on stars known to or suspected of harboring protoplanetary disks as well as some debris disks. Here 
we present some results illustrating the performance of the instrument. Actual science results and analysis will 
be the topic of forthcoming papers. The direction of the scattering polarization from a disk is perpendicular to 
the direction to the star. In the uncalibrated double-difference images, this polarization signature produces a 
double-lobed butterfly pattern, illustrated in Figure 6. 




Figure 6. Uncalibrated fractional polarization images of the AB Aurigae protoplanetary disk illustrating the butterfly 
signature produced by the tangential polarization pattern. For comparison, the image on the right shows the diskless star 
HD112815. The images are clipped at ±1% degree of polarization. 

To demonstrate the advantage of the beam-exchange method, we have generated uncalibrated polarization 
images using only a subset of the data, emulating the performance of a single beam instrument and an instrument 
that does not include modulation. The results are shown in Figure 7. 

In the left image, where the left and right beams are combined without modulation, differential aberrations 
between the two beams clearly destroy the butterfly pattern. In the middle image, in which subsequent modulated 
images from a single beam are used, the change in atmospheric seeing between the subsequent images produce 
smearing and a less defined butterfly pattern. The sharpest polarization pattern is clearly produced by the 
leftmost image, produced using the beam-exchange method. 

The ExPo polarimeter is not intended primarily to provide extremely accurate absolute polarization measure- 
ments, but to be very sensitive to spatial variations in the polarization. Still, the calibration procedure described 
in the previous section yields results accurate down to the 1% range. Table 4 presents typical accuracies when 
back-calibrating the dome flats taken at different calibration polarizer angles: 

Observations of diskless and unpolarized standard stars must be performed in order to check that the in- 
strumental polarization introduced by the Nasmyth flat is correctly compensated. Table 5 shows the calibrated 
Stokes vectors of two diskless stars, demonstrating that this is the case. 
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Figure 7. The advantage of the beam-exchange method. Left image: dual-beam without modulation. Middle: Single 
beam with modulation. Right: dual beam and modulation. Target is AB Aurigae. 

Table 4. Typical dome flat back-calibration residual calibration errors. 



Calibrator 


Using FLC & 22.5 


position 


Error in Q 


Error in U 


+Q 


0.17% 


1.03% 


+u 


0.22% 


0.20% 


-Q 


0.16% 


1.00% 


-u 


0.22% 


0.20% 


Using FLC 45 & 67.5 


+Q 


0.24% 


0.63% 


+U 


0.39% 


0.34% 


-Q 


0.22% 


0.59% 


-u 


0.39% 


0.34% 



A more visual confirmation can be obtained by simply mapping the polarization of tlie central part of a 
diskless star. The result can be found in Figure 8. 

A final analysis that can be performed using the diskless stars is to check how well the instrument rejects 

unpolarizcd light as a function of the radial distance to the star. Figures 9 and 10 show a radial profile of the 
achieved contrast and fractional polarization based on azimuthal averages. Figure 9 shows the total contrast: the 

Table 5. Calibrated linear Stokes vector elements for two diskless stars. The values are averages of a 5 x 5 pixel area 
centered on the star. 



Calibration 


Using FLC & 22.5 


star 


Measured Q Measured U 


HD86275 


-0.015% 0.002% 


HD107170 


-0.024% -0.023% 


Using FLC 45 & 67.5 


HD86275 


0.047% -0.002% 


HD107170 


0.014% 0.017% 



10 



HD122815 




-6 



-4 



-2 







2 



4 



6 



[arcsec] 



Figure 8. Polarized intensity image of the unpolarized reference star HD122815 overlaid with calibrated polarization 
vectors. The reason we still see a signal is that the polarized intensity, defined as \/ + U^, includes the absolute value 
of the photon noise. The vectors show no preferred orientation or other structure, showing that there is no residual 
instrumental polarization after the calibration. 

polarized flux divided by the peak total stellar flux. The contrast achieved is 10~* at the star position, increasing 
to 10^^ at about 4" from the star. It should be noted that fewer frames were taken in the observation of this 
diskless star in comparison with regular science targets. The higher photon noise limits the contrast achieved. 
This is clarified by the Figure 10 showing the fractional polarization (polarized flux divided by total flux at 
the present radial distance from the star) plotted together with the profile of the photon noise detection limit. 
The measured fractional polarization closely traces this theoretical detection limit, except very close to the star, 
where the statistics become degenerate because of the small number of pixels that can be azimuthally averaged, 
and far away from the star, where there is hardly any signal detected and the result becomes very sensitive to 
small flat-field or dark correction errors. The fractional polarization that can be detected is therefore limited by 
the photon noise. Higher contrasts can be achieved using longer exposure times. 

Having shown that our sensitivity is limited by the photon noise, it is clear that high-contrast polarimetry is 
a photon-hungry application. So while the polarization sensitivity and accuracy are key performance parameters 
for a polarimeter, another important characteristic is the efficiency with which these results are achieved, i.e. 
how well the instrument makes use of the available photons. [19] define the polarimetric efficiency with which a 
polarimeter measures the Stokes elements as: 



With Y being the inverse of the signal matrix X as described in the previous chapter, m the number of 
measurements and i the index of the Stokes element. This expression makes no assumptions on how many 
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(7) 
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MWC480 Tangential Contrast tanP/max(l) 




Radiol distance (orcsec) 

Figure 9. Contrast reached for increasing distance from the star. At the peak of the star, the contrast reached purely by 
polarimetry is 10~^. This increases further out as the stellar intensity decreases. This target was observed for about 10 
minutes, higher contrasts can be reached with longer observations. 

MWC480: fractional polarization vs photon noise 



1 2 3 4 5 6 

Radial distance (orcsec) 

Figure 10. Fractional polarization averaged in circles around the star for increasing radial distance (solid line). The 
sensitivity is limited by photon noise and thus reaches its maximum close to the peak of the stellar flux. Dashed line: 
photon noise limit per pixel, averaged in circles around the star. 

measurements are performed to determine the Stokes element and how these measurements are performed. It is 
thus a global measure of the polarimeter efficiency. The efficiencies per Stokes element are defined so that: 



4 + 4 + 4 < 1 



(8) 



As we do not measure circular polarization, the efficiency in V can be ignored. If we then require equal 
efficiency in Stokes Q and U we see that the maximum theoretical efficiency that can be reached is -\/(L5 = 0.707. 
Table |6] presents the efficiencies calculated for the individual beam and FLC state combinations (not the double 
difference) normalized with this theoretical maximum. For the individual beam / FLC state combinations, all 
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four FLC angles are used to generate the calibration matrices. 



Table 6. ExPo polarimetric efSciency per stokes element for individual FLC state & beam combinations. 



FLC state 


Beam 


Q 


U 


A 


Left 


0.79 


0.79 




Right 


0.81 


0.82 


B 


Left 


0.79 


0.79 




Right 


0.79 


0.81 



6. SAMPLE RESULT: AB AURIGAE 

As a sample science result, we present our imaging polarimetry observations of the well-known circumstellar 
disk around AB Aurigae. Imaging polarimetry is an excellent tool to observe the faint, linearly polarized light 
scattered by the circumstellar disk material while suppressing the unpolarized starlight. [20] and [21] have presented 
early studies of circumstellar disks using this technique. 



The protoplanetary disk around AB Aurigae, a 4±1 Myr old Herbig star, has been shown 22 23 to have 



several spiral arms. Imaging polarimetry observations by 24 revealed an annulus with a gap at about 100 AU 
from the star. This gap was speculated to be caused by a companion, but 25 showed that the structure can 
also be explained by the geometrical projection effect of the inclined disk. Only recently 26 presented very high 
resolution imaging polarimetry of the AB Aurigae disk. Their results reveal detailed structure in the inner region 
(up to 150 AU) of the disk, as well as clearly revealing the spiral arms extending out to about 450 AU. These 



results were all achieved using AO-assisted (or space-based in the case of 25 ) coronagraphy in the near-infrared. 

Our observations date from January 2"^^ 2010. The imaging has been performed without the use of a color 
filter, thus spanning the full 400-900 nm wavelength range of our detector. Other observations have been 
performed using a Schott OG590 filter, but the results are inferior. This is not surprising given the AOe spectral 
type of AB Aurigae and the higher scattering efficiency towards the blue end of the spectrum. The average 
seeing during the observations was 1.0". 

Figure 11 shows two calibrated polarized intensity images, each generated by combining two sets of observa- 
tions with the position of the FLC modulator 22.5° apart. The structure in both images is very similar, even 
thought the seeing was slightly worse for the observations used in the right image. The disk is slightly elongated 
along the NE-SW axis, which is consistent with previous observations. In the inner region, up to about 1" from 
the star, we detect and elongated, almost bipolar structure oriented in the perpendicular direction. This struc- 
ture is about at the highest resolution that can be resolved in our seeing-limited observations, but corresponds 



well to regions with higher polarized intensity to the NW and SE in the observations by 26 . Intriguingly, the 
first black contour in the left image appears to trace a similar annulus with a gap as reported by f24', although 
the position is rotated slightly. Again, this is at the resolution limit and does not appear in the right-hand image, 
but both show traces of an annulus with a cleared inner region to the north of the star. 

What can be seen very clearly in both images is a spiral arm starting at about 150 AU ( 1") east of AB Aurigae 
and curving outward in the clockwise direction. To bring this structure out more clearly we have azimuthally 
averaged the image and then subtracted this from the original, producing Figure 12. The arm labeled 'SI' by 



26 shows up very clearly, while there are hints of two other arms. 

Figure 13 shows the disk polarized intensity overlaid with a vector plot mapping the polarization orientation. 
This is a useful test to ascertain that it is indeed the scattering polarization that is observed and that there are 
no offsets caused by errors in the calibration or the removal of instrumental and sky polarization. As expected, 
the polarization direction is perpendicular to the scattering plane, producing the azimuthal vector plot. 
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Figure 11. Broadband-visible polarized intensity images of AB Aurigae. Observations at two positions of the liquid 
crystal polarization modulator have been combined to create each calibrated image: 0° and 22.5° for the left image, 45° 
and 67.5° for the right image. The color scale is logarithmic and histogram-equalized. Black contours go from 80% of the 
maximum polarized intensity downwards in steps of 20%, white contours go from 10% of the maximum down in steps of 
2%. 

7. EXPO UPGRADES 

An adaptive optics (AO) system has been developed for ExPo over the last year. This was done for a number of 
reasons: First of all, as the polarimetry is performed using large numbers of short-exposure frames, the increased 
image stability is expected to increase the polarization sensitivity significantly. The increased stability also means 
that the coronagraph built in to ExPo can now be used more effectively. Finally, the AO will hopefully increase 
the spatial resolution, and thus the level of detail we can see in the circumstellar structures. Simulations with 
the AO package PAOLA have shown that for 1" seeing and a 12*'' magnitude star, an AO with a 50-actuator 
DM would improve the Strehl of the target by 20-25 times. 

We have developed a custom design, intended to affect the polarimetry as little as possible. The adaptive 
optics system uses a 97-actuator deformable mirror (DM) from ALPAO. Two off-axis parabolas at reflection 
angles of 30° are used to form a pupil image on the DM and then form a secondary focus. Avoiding 45°- 
reflections significantly decreases the instrumental polarization introduced. A non-polarizing beamsplitter is 
used to split the light between the science arm of the instrument and the Shack-Hartmann wavefront sensor. 
Figure [14] provides an impression of the optical and mechanical setup of the AO system. 

The ExPo AO system had its successful commissioning run in the summer of 2012. 

A second upgrade currently being worked on is a spectro-polarimctric integral field unit for the ExPo instru- 
ment. This is a technology demonstration project for the future EPICS-EPOL cxoplanet imaging polarimeter 
for the E-ELT. The IFU will use a lenslet-based design to sample the image plane and generate low-resolution 
spectra for each sub-image. This will combine the high contrast provided by polarimetry with the ability to 
obtain spectral information of various points/areas of the target observed. 

8. DISCUSSION & CONCLUSION 

The results we have presented demonstrate that our choice for the fast-modulating beam-exchange method 
and our relatively simple polarization-friendly optical design results in a sensitive imaging polarimeter, able 
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Figure 12. At left, the azimuthal average has been subtracted from the AB Aurigae polarized intensity image. At right, 
the image has been smoothed with a box-car width of 2 pixels. At least one of the spiral arms can clearly be identified. 
The arms have been labeled with the corresponding numbers from [26]. The color scale is logarithmic. 
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Figure 13. Polarized intensity image of AB Aurigae overlaid with a vector plot of the polarization direction clearly 
showing its azimuthal structure. 

to produce high-contrast observations of circumstellar environments without the aid of an AO system or a 
coronograph. While the sensitivity of around 10^^ is not as high as some aperture-integrated polarimeters such 
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Figure 14. CAD-drawing of the ExPo Adaptive Optics system. 



as PlanetPol 27 and POLISH 28 , the spatial information is highly useful, not only because of its inherent 
potential to detect and study structures in circumstellar envelopes but also as a means to check the observations 
for global errors introduced by instrumental or sky polarization. The calibrated polarization accuracy of 1% that 
can be reached with ExPo is similar to these instruments. 

Being seeing-limited, the spatial resolution achieved with ExPo is not anywhere near that which can be 



achieved by current (e.g. HiCIAO: [29]) and future (e.g. GPI: [30], SPHERE-ZIMPOL: [31]) AO-assisted 
coronographic imagers on 8-m telescopes. But none of these instruments have been designed specifically with 
polarimetry in mind and are likely to suffer sensitivity and accuracy losses incurred by the many reflections before 
the light arrives at the polarimeter module. Moreover, with the exception of Sphere-ZIMPOL, these instruments 
operate in the near-infrared, forgoing the inherently higher resolution and higher scattering efficiency available 
at shorter wavelengths. 

With our sensitivity limited by photon noise, ExPo would of course benefit from operating at a bigger telescope 
or adding up more frames for a longer accumulated exposure time. The last option is however complicated by 
seeing variations that occur over the longer timescale. Apart from the obvious attraction of higher resolution, 
this is another motivation for developing a polarization-optimized AO system for ExPo. An AO system will also 
allow us to make better use of the coronograph included in ExPo. 

The suppression of starlight achieved by combining coronography and imaging polarimetry, coupled to the 
resolution and photon-flux provided by ELT-class telescopes and their associated extreme-AO systems, will 
enable contrasts of 10^^ - 10"^ to be achieved, the threshold at which Earth-like planets will be come detectable. 
We are currently applying the experience gained with ExPo to the design of the EPOL imaging polarimeter arm 
of the EPICS planet imaging & characterization instrument for the European ELT. 
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